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The stability of graphene nanoribbons in the presence of typical atmospheric molecules is system- 
atically investigated by means of density functional theory. We calculate the edge formation free 
energy of five different edge configurations passivated by H, H2, O, O2, N2, CO, CO2, and H2O, 
respectively. In addition to the well known hydrogen passivated armchair and zig-zag edges, we 
find the edges saturated by oxygen atoms to be particularly stable under atmospheric conditions. 
Saturation of the zigzag edge by oxygen leads to the formation of metallic states strictly localized on 
the oxygen atoms. Finally, the vibrational spectrum of the hydrogen and oxygen passivated ribbons 
are calculated and compared. 

PACS numbers: 73.20.At, 73.22.Pr, 71. 15. Mb 



I. INTRODUCTION 

Graphene, a single atomic carbon layer, has attracted 
extensive research interests in the last few yearsi - — . Re- 
cent advances in experimental techniques have lead to 
the fabrication of several graphene-based structures for 
both fundamental physics investigations and promising 
applications. Graphene nanoribbons (GNRs), in partic- 
ular, have received considerable attention due to their 
possible role in future carbon-based nanoelectronics^. 
GNRs are graphene strips of finite width, in which lat- 
eral confinement opens an electronic gap as opposed to 
the vanishing gap in infinite graphene sheets. 

Theoretical calculations and experiments have shown 
that the electronic and transport properties of GNRs are 
strongly influenced by the actual atomic configuration of 
the edges 7 which, contrary to bulk graphene, are very 
reactive. A detailed knowledge of the stability and elec- 
tronic properties of the different types of edges is there- 
fore required in order to understand and ultimately de- 
sign GNRs with specific properties. 

Most theoretical calculations on GNRs have consid- 
ered the the armchair and zigzag edges with or without 
hydrogen passivation. However, recent calculations have 
shown that other edge configurations can also show high 
stability 8 . Since molecular hydrogen is one of the most 
common contaminants even in ultra high vacuum exper- 
iments, it is certainly relevant to consider its influence 
at the edges 9 . On the other hand, especially at room 
temperature and under ambient conditions, the effect of 
other gas molecules on the edge stability should also be 
considered. This is also relevant from the perspective of 
chemical functionalizatio n 10 ! 11 . 

In this paper, we study the stability of five different 
edge geometries in the presence of H2, O2, N2, CO, CO2 
and H2O using first principles calculations. The different 
edge geometries for the pristine ribbons with the corre- 
sponding abbreviations are shown in Fig. [TJ We find that 
the edges passivated by atomic oxygen are the most sta- 
ble under atmospheric conditions closely followed by the 



zig-zag and armchair edges with hydrogen adsorbed. The 
oxygen passivated zig-zag GNR is found to be metallic 
with the states at the Fermi level strictly localized on the 
oxygen p orbit als. The vibrational spectrum of the hydro- 
gen and oxygen passivated GNRs are calculated and dis- 
cussed. Throughout the paper we focus on free-standing 
graphene and thus neglect any effect related to the sub- 
strate. Whether this is a reasonable assumption clearly 
depends on the strength of the substrate-graphene inter- 
action. This is likely to make the edges less reactive to- 
wards gas species, since the edges would be already stabi- 
lized by the chemical bond to the substrate. For metallic 
substrates the nature of the substrate-graphene interac- 
tion is still unclear, but DFT calculations indicate that 
it can vary considerably depending on the meta l 12 i 13 . 



II. COMPUTATIONAL DETAILS 



Density functional theory (DFT) calculations are per- 
formed with the GPAW cod o 14 i 15 which is a real space 
implementation of the the projector augmented wave 
method^. The widths of the nanoribbons are in the 
range 20 — 23 A, corresponding to 10-12 carbon dimers, 
which is needed in order to get fully converged results 
with respect to the width. Note that this size is also 
relevant to modern experiments, where sub-lOnm wide 
GNRs are now achieved^. The GNRs are separated 
by 10 A in all directions. All calculations are spin- 
polarized and use the RPBE 17 exchange-correlation func- 
tional. A grid-spacing of 0.18 A is used, and a (1x1x10) 
Monkhorst-Pack /c-point grid is employed to sample the 
Brillouin zone along the periodic direction of the nanorib- 
bons. We relax all the structures until the maximum 
force is lower than 0.05 eV/A. 
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FIG. 1: The standard and reconstructed edge geometries (and 
corresponding abbreviations) considered in this study. The 
figures show two unit cells of the nanoribbons with periodicity 
along the z-direction. 
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TABLE I: Numerical values of the adsorption energies, -E a ds, 
corresponding to Fig. [2] All values are in eV. 



III. STABILITY 

A. Adsorption energies 

For each of the five edge configurations we calculate the 
adsorption energy of the different molecules according to 
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where ^; ribb + mo1 i s the total energy of the ribbon with the 
adsorbed gas molecule, iV m oi is the number of adsorbed 
molecules, E mo1 is the total energy of the molecule in the 
gas phase and £; ribb is the total energy of the pristine 
ribbon. In order to avoid the erroneous DFT description 
of the O2 gas-phase triplet groundstate, we use H2O and 
H2 as suggested in Ref. [HI. The adsorption energies are 
summarized in Fig. [2] and shown in HI 

It can be seen that saturation by oxygen is the most 
favorable for all edge configurations. We note that this 
is in agreement to Ref. 19 in which the authors conclude 



FIG. 2: (Color online) Adsorption energy, i? a ds, for the dif- 
ferent molecular species on the different edge configurations 
in eV. In the case of nitrogen and oxygen we have considered 
both the molecular and dissociated forms on the edge. In all 
cases the adsorption energy is given per molecule. 



that zig-zag nanoribbons are more likely to be oxidized 
than hydrogenated. 

While the high reactivity of the zig-zag edge is well 
estabilished 20 , we also find the ac(56) to be very reactive; 
in fact for several species this is the most reactive edge 
configuration. As for the zig-zag edge, the origin of the 
high reactivity of the ac(56) edge is a (spin polarized) 
peak in the density of states close to the Fermi energy. 
Consistent with their high reactivities, the pristine zz 
and ac(56) edges also have the highest edge formation 
energies (see Table HI]) , i.e. they are the least stable edges 
in the absence of molecular gas species. 

Finally, we note that all the reported results are for 
adsorption at the edge carbon atoms since these are by 
far the most reactive sites. For example we find that H 
adsorption on a carbon atom located next to the edge is 
already unfavorable compared to H2 in the gas phase. 



B. Edge formation energies 

The adsorption energies discussed in the previous sec- 
tion are relevant once a given edge configuration has been 
formed. More generally, to thermodynamic stability of an 
edge in the presence of gas molecules should also take the 
energetic cost of forming the edge into account. To cal- 
culate the edge formation free-energy we first calculate 
the formation energy of the pristine edges, 
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where E rihh is the total energy of a ribbon with Nq car- 
bon atoms in the supercell and £; bulk is the total en- 
ergy per atom in bulk graphene. Our calculated energies 
shown in Table [II] are in very good agreement with Ref. 
d. We remark, that for pristine GNRs the zz(57) recon- 
struction has the lowest edge energy. 

The adsorption energy E a( ±s is used to obtain a free- 
energy by taking into account the chemical potential, 
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TABLE II: Edge formation free-energy at OK and 300K for 
the different nanoribbons and adsorbed species. All energies 
are in eV/A. 



FIG. 3: (Color online) Edge formation free energy of graphene 
nanoribbons with different edge structures (different symbols) 
and different adsorbed molecules at the edge. The free energy 
is evaluated at 300i^. The pressure for each structure with 
a given adsorbate is taken to be the partial pressure of the 
gas-phase adsorbate at atmospheric conditions. 



Mgas^ of the gas at a given temperature and pressure 

G ads = -£ ads - PgasMgas (3) 

where p gas = N mo \/L is the density of adsorbed species. 
The chemical potential at temperature T and partial 
pressure P is given by 

= H°(T) - H°(0) - TS°(T) + k B Tln (J^J (4) 

where H° and S° are enthalpies and entropies of the gas 
phase molecules at P° = 1 bar. 

Finally, the edge formation free-energy in the presence 
of adsorbates is obtained as 
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The vibrational contributions to Gf, a ds are not included 
since they are expected to be insignificant 21 . 

Fig. [3] summarizes the calculated edge formation free- 
energies at 300K and atmospheric partial pressures for 
the different gas species. The qualitative picture does 
not change if the temperature and pressure ranges are 
varied within experimentally accessible ranges. The re- 
sults are also reported in Table HI1 together with the edge 
formation energies calculated at OK. It is immediately 
apparent that edges saturated by atomic oxygen are the 
most stable. In particular, the lowest energy configura- 
tions are the standard zigzag and armchair, followed by 
the reconstructed zz(57); the other reconstructed edges 
are generally less stable. This agrees with Ref. 9 where 
hydrogen passivation was considered. The armchair and 



zigzag edges saturated by hydrogen are also quite sta- 
ble thus making H and O saturations the most stable 
ones. For this reason we will focus in particular on the 
oxygen saturated GNR, i.e. the most stable one, in the 
rest of the paper. Note that the negative formation free- 
energies for some of the oxygen saturated edges mean 
that a graphene sheet will spontaneously break in or- 
der to generate edges and lower its energy. However one 
needs to be careful in analyzing these results, since signif- 
icant energy barriers may be involved in the process and 
they are not included in our study. Furthermore, once 
the oxygen molecule has dissociated, atomic oxygen is 
also reactive on the planar graphene sheet (3.23 eV com- 
pared to atomic oxygen 22 ), thus making the dynamics 
of the dissociated oxygen potentially very complicated. 
Oxygen saturated edges were also found to be the most 
stable in another study 23 , where the authors also con- 
sider edge configurations beyond saturation by a single 
atom. 

In addition to H and O, the adsorption of H2O and CO 
is also found to stabilize most edges. In the case of H2O, 
the dissociation of the water molecule is always energeti- 
cally favoured compared to the molecular adsorption. In 
particular, it dissociates into OH and H on ac, ac(56), 
zz(57) and ac(677). On the zz edge, instead, H2O spon- 
taneously dissociates into an adsorbed oxygen at the edge 
and releasing the two hydrogen atoms in the gas phase. 
Nitrogen species and CO2 bind rather weakly or not at 
all on all edge configurations. 

These calculated data can furthermore be used in order 
to construct phase diagrams as a function of the chemical 
potential of a given gas species. This might be relevant 
for guiding experimental investigations in which GNRs 
are studied in the presence of a given gas species. An 
example is shown in Fig. 2] for the case of dissociated 
oxygen at the edges. This shows that for any experi- 
mentally accessible pressure of O2, the most stable edge 
configuration is the standard zig-zag, saturated by oxy- 
gen. 
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FIG. 4: (Color online) Edge formation free-energy for the 
different edge configurations of a graphene nanoribbon in a 
background of O2 gas. 




Wave vector 

FIG. 5: (Color online) Band structure of a 22 A wide zig- 
zag GNR saturated with O. The two insets show the states 
corresponding to the bands crossing the Fermi level, marked 
with black circles. 



IV. OXYGEN SATURATED ZIG-ZAG EDGE 

In the previous section we found that the zig-zag GNR 
saturated by atomic oxygen is the most stable thermody- 
namically. For this reason, we investigate in this section 
the ground state and vibrational properties of the oxy- 
gen passivated zig-zag GNR and compare it to the well 
studied hydrogen passivated case. 



tional. Contrary to the hydrogen saturated case, its 
ground state is spin-paired and metallic. The states at 
the Fermi level are localized at the edges, and they have 
weight on the p orbitals of the oxygen atom. One state 
has p z symmetry and it decays into the GNR away from 
the edge with weight on only one carbon sublattice (left 
inset of Fig. [5]). The other state is almost completely 
localized on the out of plane p oxygen orbital (right inset 
of Fig. [5]). This is also in constrast to the hydrogen sat- 
urated case, where the edge states are entirely localized 
on the edge carbon atoms. 



B. Vibrational properties 

Low energy spectroscopies probing the vibrational ex- 
citations have traditionally been a powerful tool for struc- 
ture determination in molecules. It has recently become 
possible to study differences in the edge orientation of 
GNRs using Raman spectroscop y 24 ^ 5 . In view of this we 
have calculated the vibrational spectrum of the zig-zag 
GNRs passivated by hydrogen and oxygen in order to 
identify features specific to these types of edges. 



1. Computational details 

To obtain the phonon frequencies we first compute the 
dynamical matrix in real space. The latter is obtained 
by numerical differentiation of the forces. A careful pro- 
cedure to obtain the equilibrium configuration is found 
essential in order to obtain an accurate description of 
the phonon frequencies. Thus, for a single unit of the 
zzGNR, we used a grid spacing of 0.1 A and a /c-point 
sampling of (1 x 1 x 10). The relaxation of the equilibrium 
structure is continued until the maximal residual forces 
are less than 0.02eV/A. 

The dynamical matrix is evaluated in a supercell con- 
taining seven primitive cells of the GNR. All atoms in one 
of the primitive cells are displaced by ±0.02 A in the x, y, 
and z directions (a total of 6 x N calculations ) and the 
forces on all remaining atoms within a truncation radius 
of 7 A of the displaced atom are stored. This requires 
a total of 6 x N calculations where N is the number of 
atoms in a primitive cell. For these calculations we use 
a grid spacing of 0.18 A and a 1 x 1 x 3 /c-point sam- 
pling. We have verified that our results have converged 
with respect to these parameters and the truncation ra- 
dius. The real space dynamical matrix is obtained from 
a finite difference of the forces, and the eigenvectors and 
frequencies are obtained after Fourier transformation. 



Bandstructure 



Fig. [5] shows the Kohn-Sham band structure for the 
oxygen saturated GNR calculated with the RPBE func- 



2. Results 

The GNRs are finite in their width direction and thus 
they are subject to boundary conditions similar to that 
of a string with free ends. The edges therefore only allow 
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standing waves perpendicular to the ribbon axis and the 
phonon wave vector for a GNR consisting of N carbon 
dimers will be quantized as q± y7l = irn/w, where w is the 
width of the dimers and n = 0,...,iV — 1 each correspond- 
ing to a normal mode. The dispersion relations there- 
fore consist of six fundamental modes each with (TV — 1) 
overtones. A study of the T-point spectrum according to 
this classification has been done in a previous study^ for 
hydrogen saturated nanoribbons with purely zigzag and 
armchair edges. The modes can be separated into the 
three categories according to their longitudinal, trans- 
verse, and out-of-plane nature. It is straightforward to 
identify which category a particular mode belongs to by 
weighting each of the three different spatial entries in the 
polarisation vector relative to its total norm. Each of the 
spatial categories can be separated into acoustical and 
optical parts giving rise to six series. The T-point spec- 
trum is found to be similar to hydrogen saturated ribbons 
where a splitting into optical ( 160-200) me V and acousti- 
cal modes (0-160)meV is observed, as shown in Fig. [6j 
The transverse acoustical modes are seen to follow a near 
perfect 1 / t/w behaviour which is understood in terms of 
the ribbon acting as a spring perpendicular to its axis 
with fixed spring constant. 
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FIG. 6: The longitudinal (left) and transverse (right) r-point 
spectrum of the oxygen saturated zigzag nanoribbon as func- 
tion of the number of carbon dimers. The fundamental mode 
of the acoustical series is zero for all widths and is not shown. 

A separation of modes associated with the saturation 
atoms are not as simple in the case of oxygen as it is for 
hydrogen. Contrary to the hydrogen passivated ribbons 
where three degenerate pairs of edge localized modes are 
found at 380meV, 150meV and 125meV, only one local- 



ized degenerate pair of carbon-oxygen modes shows up, 
namely the degenerated stretching modes (transverse di- 
rection). Similarly to the hydrogen passivated case, this 
pair is found to be independent of the ribbon width and 
is located around 182meV, except for the very narrow 
ribbons (N < 4) where it is slightly higher. The last 
four modes are not separated from the remaining modes 
as for the case of hydrogen. In fact, the oxygen atoms 
participate as a continuation of the carbon system and 
enters in the bulk vibrational modes giving rise to two ad- 
ditional acoustical modes in the longitudinal series and 
the out-of-plane series. This coupling to the acoustical 
modes is understood in terms of the comparable masses 
of the two species. Also, for the transverse series one 
finds that an optical mode has been converted into an 
additional acoustical one. The optical modes of the lon- 
gitudinal and out-of-plane series on the contrary remains 
similar to what is found for the hydrogen saturated rib- 
bons, that is, the series are constituted of normal modes 
only coming from the carbon atoms. 
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FIG. 7: (Color online) The total vibrational density of 
state of the 16 dimer wide (33A) hydrogen(red) and oxy- 
gen (blue) saturated zigzag nanoribbon in comparison with 
graphene (black). The density of states weighted on the sat- 
urating atoms shown for the hydrogen and oxygen saturated 
16 dimer zigzag ribbon, respectively. 

The total vibrational density of states of the 16 dimer 
wide zigzag ribbon is obtained by integration in the 
TKM-direction and is shown in Fig. [71 Note that it 
represents the whole spectrum of states and not just the 
T-phonons. Except for some distinct peaks the total den- 
sity of states for the hydrogen and oxygen saturated rib- 
bons represents that of graphene fairly well. In order to 
compare the deviations in the density of states associ- 
ated with the edge, the states localised on the edges have 
been separated out for each of the saturating species. 
The peaks present in the edge contribution are to a large 
degree recognisable in the corresponding total density 
of states spectrum. For the hydrogen saturated rib- 
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bon, the two peaks observed are the vibrational mode 
perpendicular to the carbon-hydrogen bond and out-of- 
plane mode. For the oxygen saturated ribbon only the 
transverse carbon-oxygen has a truly edge-localized vi- 
brational pattern while the peaks observed in the acous- 
tical part of the spectrum originates from individual nor- 
mal modes in the longitudinal and out-of-plane series. 

V. CONCLUSION 

As a general reminder, it should be noted that the 
results in this study do not consider the effects of the 
substrates. In most experiments, GNRs are grown or 
deposited on a substrate, which might affect the energet- 
ics at the edges. The interaction of the substrate with 
graphene edges is likely to be much stronger than the 
interaction with the bulk graphen o 12 i 13 . 

In conclusion we have studied the adsorption of com- 
mon gas molecules at the edge of several types of GNRs 
using density functional theory. We have calculated the 
edge formation free energy in the presence of the adsor- 
bates at the different edges and shown that passivation 
by atomic oxygen is the most stable closely followed by 
the hydrogen passivated GNRs. In contrast to the hydro- 
gen saturated zig-zag edge, passivation by oxygen leads 



to a non-magnetic and metallic groundstate. 

For the vibrational properties, the comparable masses 
of oxygen and carbon lead to a larger degree of mixing 
between localized edge and bulk modes. In particular, it 
is found that the oxygen atoms act as a continuation of 
the carbon system and participates actively in the lon- 
gitudinal and out-of-plane series. Only the stretching 
(transverse) mode of the carbon-oxygen bond is found to 
be completely localized at the edge and separated from 
the normal mode series. 

Since O2 and H2 are almost always present under re- 
alistic conditions, our results suggest that it might be 
very difficult to functionalize the edge atoms with dif- 
ferent species. We suggest that the substitutional route, 
for example via electrothermal reactions as demonstrated 
in 27 , might turn out to be an easier approach for doping 
GNRs. 
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